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Abstract 

We investigate the quantum information reclaim from the environment after amplitude damping has 
occurred. In particular we address the question of optimal measurement on the environment to perform 
the best possible correction on two and three dimensional quantum systems. Depending on the dimen- 
sion we show that the entanglement fidelity (the measure quantifying the correction performance) is or 
is not the same for all possible measurements and uncover the optimal measurement leading to the max- 
imum entanglement fidelity. 

PACS number: 03.67.Pp, 03.65.Yz 



1 Introduction 

Decoherence is a fundamental problem for the realization of quantum information tasks. The unavoidable 
interaction of a system with the surrounding environment disturbs the state and causes decoherence [1]. This 
situation is described by a unitary interaction between the system and the environment. The state change of 
the system can be described by a CPT map (completely positive trace preserving map). Usually this implies 
an irreversible flow of information from the system to the environment with a consequent wash out of coher- 
ence features of the system. Having access to the environment one could think to reclaim the information 
lost into the environment and restore the (quantum) coherence features of the system with a proper action 
on it. This would be a typical example of feedback (closed loop) control [2] in which the actuation on the 
system would be based on classical information gathered from environment. This method is practical for 
systems interacting with sufficiently controllable environments, that is, environments on which the neces- 
sary measurements can be performed (e.g. in experiments inside cavities). 

Recently, a lot of attention has been devoted to the feedback control scheme from different aspects. In 
[3, 4] the capacity for this scenario has been studied and in [5] it has been shown that in the limit of infinite 
number of iteration of a map on qubits and qutrits, the effect of decoherence can completely be removed. In 
[2] it has been shown that if the initial state of the environment is pure, then any measurement on the envi- 
ronment corresponds to a specific Kraus representation of the map describing the evolution of the system. 
For a given measurement the best recovery scheme has also been introduced [2]. However, it is clear that 
the success of this technique resides on the optimality of measurement. 

Here we address the question of what kind of measurement can provide the highest achievable coherence 
after performing the optimal recovery. In particular we consider the nontrivial case of amplitude damping 
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channel describing the loss of energy (hence information) from system to environment. Due to the corre- 
spondence between measurement and Kraus representation of the map, any classical result of measurement 
a can be interpreted as the error (Kraus operator) occurred on the system. Since the canonical Kraus 
operators are independent of each other, a naive guess about the optimal measurement might be the one 
corresponding to canonical Kraus representation. Actually we show that for qubits all possible measure- 
ments on the environment are equivalent in the sense of recovering the errors by feedback control. We also 
show that for amplitude damping channels acting on qutrits the situation is different and the measurement 
corresponding to canonical Kraus representation is not the optimal measurement on the environment. 

The layout of the paper is as follows. In Section 2, we briefly present the main conceptual and computational 
tools needed to restore quantum coherence by means of a quantum feedback control scheme. We introduce 
amphtude damping channel in Section 3. For d = 2, we show that all decompositions or all measurement 
on the environment are equivalent for restoring quantum coherence via feedback control. In Section 4, we 
study the amplitude damping channel for qutrits. We find that for d = 3 the canonical Kraus decomposition 
is sub-optimal. Indeed, by means of analytical methods, we construct a class of (non-canonical) Kraus de- 
compositions leading to an higher performance than the one obtained via the canonical decomposition. The 
paper concludes in Section 5 with a discussion. 



2 Quantum feedback control scheme 

In this Section we describe the general scheme for error correction using the information we gain by mea- 
suring the state of environment. The evolution of state in interaction with an environment or the action of 
any quantum channel on the input state can be described by considering the unitary evolution of the system 
together with the environment. Considering Hilbert spaces of initial and final states and and also 
the Hilbert spaces for initial and final environments /Ci and /C25 the general evolution can be described by 
a unitary operator IJ : li\ ® K\ ^ ® /C2. The final state of the system is given by tracing over the 
environment after the interaction, 

T{p)=Tt^^\V{p®g)V\ (1) 

where p and a are positive trace class operators belonging to ClTi), the space of all linear operators on 
Hilbert space Ti. The map T : C{T-Li) — C(V.2) is a CPT map which describes the channel or evolution of 
a system after interaction with environment. Our aim is to gather information about the errors occurred on 
the system by measuring the environment after interaction with the system . A general POVM measurement 
on K.2 is described by Mq, g >C(/C2) where 

J2Ma = I, Ma>0. (2) 

a 

The index a labels the classical result of measurement. To understand how this classical information can 
be used to gain information about the final state of the system, let us consider an arbitrary observable 
A G C{%2) and its expectation value, 

<A>= Tr^^ {T{p)A) = Tr^Tr^^ [U{p a)U\A /)]. (3) 

Replacing I from (2) in the above equation, we get 

<A> = J2Tr^^{Tr^jU(p^a)U\A^Ma)]} 

a 

= J2'^r^.('^a{p)A), (4) 

a 

where : C{%i) — > ^{^^2) is defined as, 

Ta{p):=Tr^AU{p(^CT)U\l(^Ma)]. (5) 
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Rewriting the expectation value of ^4 in the following way, 



<A>=Y,PaTr[^Al (6) 

Pa 



we can conclude that p„ = Tr(Taip)) is the probability of getting a as the classical result of measurement 
and the density matrix ^Ta{p) as selected state of the system after measurement. Therefore by performing 
the measurement on the environment the channel is decomposed as, 

T = J2Ta. (7) 

a 

Let us recall from [2] that the most informative measurements on the environment are those for which the 
selected output of the channel after the measurement can be described by a single Kraus operator ta '■'Hi ^ 
'H2, 

Taip) = taptl (8) 

From now on we assume that the initial state of the environment is pure. It has been proved in [2] that when 
the initial state of the environment is pure, every decomposition of the channel in form (7) can be reahzed by 
a measurement on the environment. We also assume that we can perform the most informative measurement 
on the environment . Therefore designing different POVMs on the environment is equivalent to considering 
different Kraus representations of the original channel: 

T{p) = Y.^aPti Y.^ita = i. (9) 



a 



To perform the correction we design a recovery channel Ra ■ -^(^2) ^{T^i) that depends on the mea- 
surement outcomes. The state of the system after performing the correction is ^ Ra (Tq, (p) ) with probability 
Pa and the overall channel : C{7ii) — > C(J-ii) is given by 

Tcorr = ^2 ° ^^^^ 

a 

where o means composition of the maps. The closer the Tcarr is to the id map, the more successful is the 
scheme to recover quantum information. To quantify the performance of the correction scheme, we use 
entanglement fidehty [6] as measure. For a general channel $ : jO-{'H) )~-{T~i) with Kraus operators Ak, 
the entanglement fidehty is given by 

k 

where d = dimH and 1^*) G H is a maximally entangled state. We are interested in F(Tcorr), the 
entanglement fidehty of the corrected map. Given T{p) in (9) and using the Kraus representation of the 
recovery channel Ra ■ C{T-L2) ^i'Hi): 

R.{p') = Y.^tp^f^ Y.^fr1 = i^ (12) 

/3 /3 



the entanglement fidelity of the corrected channel becomes, 

1 

d2 



F{T^r) = i M^f^^c.)?- (13) 



a,/3 



Entanglement fidelity reaches identity and quantum information can receive complete correction if and only 
if for all a, tat]y, = Tal with Tq, > and J2a ''"a = 1 t^]. In particular, when diniHi = dim'H2, channels 
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with Kraus operators proportional to unitary are completely correctable. In this case the optimum mea- 
surement on the environment is the one corresponding to the unitary decomposition of the map and the 
correction scheme in very clear because the evolution of the system is reversible. 

For those cases where quantum information can not be restored completely, the correction schemes which 
give the best information preserving are known. More precisely, using the Cauchy-Schwarz inequality, it 
has been shown [2] that for every family of recovery channels Ra : C{H2) C{T-Li), the entanglement 
fidelity F[Tcarr) is such that 

F{T^orr)<^Y.^tr\tc.\f, (14) 

a 

where \t\ = VtH is the modulus of t. Moreover the recovery scheme through which the maximum can 
be attained, is also known [2]. Therefore even if quantum information can not be restored completely, the 
correction giving the highest value of entanglement fidelity can be applied. It is important to notice that 
these results are obtained for a given measurement or equivalently for a given Kraus representation of the 
channel and leave the question of the optimal measurement on the environment unanswered. 

To find the optimal measurement on the environment, we need to maximize the entanglement fidelity over 
all possible Kraus decompositions of the channel. However, to make the problem tractable, we restrict our 
attention to Kraus representations with the same number of Kraus operators as the canonical representation. 
To be more specific, we denote the upper bound in equation (14) by F{Tcorr)'- 

F{Tcorr) = ^Y.(^r\ta\f (15) 

a 

and maximize it over all possible Kraus representations with the same number A'" of Kraus operators as the 
canonical one. We start working from the canonical representation of a given map: 

N-l 

Tip) = E ^'^P^^k (16) 

fe=0 

with Kraus operators satisfying tr{CkC\i) = CkSk,k'^ Ck G Afterwards, considering a general N 
dimensional unitary operator V, we construct a new set of Kraus operators {Bk} for the channel 

N-l 

Bk=Y. ^k,iCi, (17) 

1=0 

and maximize the entanglement fideUty given by the new set of Kraus operators over the parameters of the 
unitary matrix. 

It is important to notice that instead of considering a general unitary operator in U{N) it is sufficient to 
perform the maximization over the parameters in SU (N). This is because any operator V G U{N) can be 
written as F = e^'^U with U G SU{N) and </!> G M. Consttucting the new set of Kraus operators by means 
of V, we get: 

Af-l N-l 

Bk=Yl ^k,iCi = e^'*' Uk,iCi (18) 

1=0 1=0 

However only the absolute values of BkS play a role in the entanglement fidelity, 

\Bk\ = ^J^k = £ UliUk,vCjCi,)l (19) 
1,1' 
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Thus, the entanglement fideUty does not depend on the phase e**^ and therefore we can restrict our attention 
to Kraus representations obtained by transformations in SU{N) without loss of generality. 

In the following Sections we find the optimum measurement to perform correction on amplitude damp- 
ing channel. This is a nontrivial case because this map is not a random unitary map. Therefore, it is not 
included in the category of completely correctable maps discussed after equation (13), with trivial optimum 
measurement and correction scheme. 



3 Amplitude damping channel 

The general behavior of processes with energy dissipation to the environment is well characterized by the 

amplitude damping channel. The gradual dissipation of energy or amplitude damping channel is described 
by the interaction between the system and the environment, modeled by harmonic oscillators, through the 
unitary operator: 

where % is proportional to the coupUng constant between system and environment. The operators a and 
(resp b and b^) denote annihilation and creation operators of the system (resp environment). By truncating 
the Fock basis of bosonic mode of the system to length d, we describe a qudit system. The Kraus operators. 
Cm = {nib\U\Ob), for this d-dimensional amplitude damping channel are given by 

d-l 

Cm=Y^ V'C(n,m)(l-p)"-"^p"^|n-m)(n|, m = 0---d-l, (20) 

n=m 

where p = sin^ x is the probabihty of loosing a single quantum of energy, \n) are number states and 
C(n,m) = ^ ^ ^ is the binomial coefficient. It is straightforward checking that this is the canonical 
representation 

n 

Using the definition of the Kraus operators in (20), it is clear that CmCm is diagonal: 

d-l 

ClCm = Yl C(n,m)(l -^.)"->"»(n|. (21) 

n=m 

Since ampUtude From equations (21) and (15), it is easy to find that the highest value of the entanglement 
fidehty for the canonical representation of the amphtude damping channel is, 

d-l 

FciTcorr) = ^Y{tr\Cm\f 
m=0 

^ d-l /d-l \ 2 

= -^Y^iYl ^C{n,m)(\-v)"v^ , (22) 

m=0 \n=0 / 

where the under script c stands for canonical. For systems with two dimensional Hilbert space, the entan- 
glement fideUty for canonical Kraus representation Fc{Tcorr) is given by. 



-Fc{TcQj-j-^ — 2 . (23) 

To see whether or not there is any advantage in using the canonical decomposition, we start from the canon- 
ical representation of the map: 



Co = |0)(0|-fv'l^|l)(l|, 
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Ci = VP|0)(1|- (24) 
considering the general two dimensional special unitary operator 

we introduce a general set of Kraus operators {B^} as follows: 

For this set of ICraus operators, the entanglement fideUty is given by 

MTcorr) = IY1 i^'il^rnl) = l±^^^ . (27) 

Comparing equations (23) and (27), we conclude that all representations with two Kraus operators give the 
same entanglement fidelity as the canonical one. 

To answer the question of whether or not the same result is valid in higher dimensions, we face the problem 
of diagonalizing d-dimensional matrices ^j^-Bfc to compute the entanglement fidelity. Therefore in the next 
Section we study the amplitude damping channel in three dimension. This will give us some insights on 
possible advantages (or disadvantages) of using the canonical representation. 



4 Amplitude damping channel for qutrits 

For three dimensional Hilbert space the canonical decomposition of the chaimel is described by the following 
Kraus operators: 

Co = |0)(0| + yr ^|l)(l| + (l-p)|2)(2|, 
Ci = v^|0)(l| + V2p(l-p)|l)(2|, 

C2 = p|0)(2|. (28) 
The entanglement fidelity corresponding to the canonical decomposition is 

Fc{Tcorr) = ^[(2 - p + ^/l^f + (VP + VM^ ' P)? + P\ (29) 

As explained in Section 2, in order to find the optimal measurement we maximize the entanglement fidelity 
over all possible Kraus representations with = 3 number of Kraus operators. To construct the general 
Kraus representation from the canonical one, we can restrict ourselves to the unitary operators in SU (3) 
without loss of generality. As starting working hypothesis, we consider the following two subgroups of 
SU{3), Gi and G2, leading to equal or higher entanglement fidelity than the one given by canonical Kraus 
representation. We show that the decomposition giving the highest value of entanglement fideUty can be 
constructed by means of unitary transformations in G2. 



4.1 Equi-canonical class 

The first subgroup Gi that we study is defined as. 




Gi := {C/i = a (3 \ \ \af + = 1}. (30) 
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The general Kraus operators constructed in terms of the elements of this subgroup are given by, 

Bq = Co, 

B2 = -pCi+aC2. (31) 



To calculate the entanglement fidelity FB{Tcorr) for this class of Kraus operators, 

2 
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I 2 

FB{Tcorr) = -^Y.^tr\Bk\ 



k=0 

we need to compute the eigenvalues of the B^B^s. It turns out that the non vanishing eigenvalues of bIBi 
are given by 

A„A; = ^i±^fB!, (32) 

with 



= p'^ + 3p{l - p)\a\'^ , 6 = 2p|a|V2(l-p)- (33) 



The non vanishing eigenvalues of -62-^2 are. 



X2,X'2 = ^^^^ (34) 

with 

a' = p^ + 3pil-pm\ 5' = 2p|/3|V2(l-p). (35) 

Therefore the entanglement fidelity is given by 



± Va'2 - 6'2 



FsiTcorr) = ^[i2-p+^l-pf + a + a' + b + b'] 



^[{2-p+^/r^f + (^^+y/2pii-p)y+p\ 



(36) 



Comparing (36) with (29), we see that for all the new set of Kraus operators B^ in (31), the entanglement 
fidelity FBiTcorr) in equation (36) equals the entanglement fidelity for the canonical Kraus representation, 

FsiTccyrr^ — Fq(Tq^^^. (37) 

This means there is no advantage in using canonical representation over the class of Kraus decomposition 
in (31). 

As a side remark, we notice that other sets of Kraus operators {-B^} leading to the entanglement fidelity 
FB(Tcorr ) in (36) Can be constructed. Introduce the two S'C/ (3) -operators gi and 32 as follows 
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51 = ^ 1 j 52 = I 1 I (38) 

and consider the new set of Kraus operators {B'f^} given by 

B'k = Y.(9j)k,iBi j = l,2. (39) 
I 

From (39), we notice that the only effect of the gs on the vector B = (Bq, Bi, ^2)* is that of cychcally 
permuting the vector components. 
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The class of unitary operators giving the Kraus representations {B'j^} from the canonical ones are simply the 
following left cosets of Gi: 



i«r + i^i =1}, 

lap + = 1}, (40) 

Therefore any Kraus representation obtained from the canonical one via transformations in 

Gi\JgiGi\Jg2Gi 
leads to the same entanglement fidelity Fc{Tccyrr) in (29). 
4.2 Super-canonical class 

In the previous subSection we showed there is a large class of Kraus representations leading to the same 
entanglement fidelity that can be achieved by the canonical representation. A more interesting question is 
whether or not we can design a measurement on the environment that gives enatanglement fidelity values 
higher than the canonical one. To answer this question, we consider a new subgroup G2 of SU (3), 



G2 := {C/2 = 1 I |7|2 + \8\^ = 1}. (41) 




In terms of U2, the new set of Kraus operators becomes 

Do = -fCQ + 5C2, 
Di = Ci, 

D2 = -6Co + 'fC2. (42) 
Following the fine of reasoning presented in the previous subSection, we obtain 



with. 



g = p2 + 2(l-p)|7|2, 

h = 2{l-p)\j\\ (44) 



Similarly 



where 



A; + - /2 V k - ^/W^ 



tr\D2\={ + +Vl^\6\ (45) 



k = p^ + 2{l-p)\6\\ 

I = 2(1 -p)\Sf. (46) 
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Figure 1: (color online) Entanglement fidelity for the equi-canonical class (dash line) and the maximum 
entanglement fidelity (solid line) found in the super-canonical class of Kraus representations. Dot line 
shows entanglement fidelity without performing any correction. 



Therefore the entanglement fidelity obtained from the new set of Kraus representation becomes (Appendix 
A) 

Fd (Tcorr) = F,{T,orr) + ^"^^ 1^, (47) 

with 



^ = \lWg + h + \5\Vk + l- (2 - p). (48) 
Since Q, is strictly positive (Appendix A), it follows that 

i^D (Tcorr) > F,{Tcorr). (49) 

We conclude that all the Kraus decompositions constructed using the elements of group G2 lead to entan- 
glement fidelity values higher than the one obtained by means of the canonical representation. Using the 
same arguments of the previous subSection, we can find a larger class of unitary transformations giving rise 
to Kraus representations with higher entanglement fidelity than the canonical one. Such transformations 
belong to the following left cosets of G2, 

91G2 := {( 7 0^1 |7|' + I<5|' = 1}, 
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6 



:= { -5 7 I hr + l^r = 1}, (50) 
V 7 5/ 

where gi and 52 are given in equation (38). Therefore the operators in the set 

G2U51G2IJ52G2 (51) 

give rise to decompositions that work better than the canonical one. We can find the best decomposition in 
this class by maximizing in equation (48) over the parameters defining the transformations in class (5 1). 
It follows that the maximum is achieved for I7I = |5| = 



n^nax = V2 + 2(l-p)2 - (2 - p) . (52) 
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Thus, the maximum entanglement fideUty in this class becomes 

1 



[FD{Tcorr)] r 



9 L 



5-2p + 2pV2(l -p) + 2^2(1 -p){2-2p + p^) 



(53) 



and can be achieved by decompositions arising from the canonical one by means of the following special 
unitary transformations: 



( il. n it \ 
1 

V2 ^ V2 / 
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(54) 

Figure 1 shows the entanglement fidelity for the case that no correction has been perfomed (dot line), for 
the equi-canonical class (dash line) and the maximum entanglement fidelity (solid line) that can be attained 
in the super-canonical class. 



4.3 Maximum entanglement fidelity 

It is impractical to analytically prove that the entanglement fideUty in equation (53) is the global maximum 
not just the local maximum in the super-canonical class. However, we are able to show this numerically. 
We generate n = 10^ sets of random unitary operators using the fact that an arbitrary operator in SU (3) is 
generated by Gell-Mann matrices, {Aj} [7]: 

8 

U = exp[— i ^ «i Aj] , (55) 

where a, are real coefficients. The chosen cardinaUty n of the randomly generated distinct Kraus representa- 
tions is the smallest positive integer number necessary to obtain convergence to the numerically found global 
maximum. Any other randomly generated set of Kraus representations with higher cardinality m, > n con- 
verges to the same global maximum. Our numerical analysis implies that such global maximum coincides 
with the analytical expression in equation (53). Thus, we conclude that [Fj:,{Tcorr)\max in (53) is indeed the 
global maximum. 



5 Discussion 

The main emphasis of this paper is to reclaim quantum information lost to the environment surrounding the 
system. Although the evolution of the system is irreversible, but having access to the environment enables 
us to restore quantum features of the system. The main idea is to perform appropriate correction on the 
system regarding the classical results obtained by measuring the environment.This scheme is an example of 
feed back control. Here we have addressed the important question of what kind of measurement provides us 
with the highest recovery performance. 

To answer this question, we have applied quantum feed back control scheme to amplitude damping channel 
describing the loss of information to the environment. After describing the amplitude damping channel in 
arbitrary dimension we have concentrated on two dimensional systems. We showed that all measurements 
corresponding to Kraus representations of the map with two Kraus operators (measurements with two out- 
comes) give the same value for the entanglement fidelity. It implies that for the amphtude damping channel 
in d = 2, the reduction of decoherence using the quantum feedback control scheme does not depend on the 
details of the measurement made on the environment. Our numeric studies show that the same statement is 
valid even if we consider measurements with three or four outcomes. 
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Although in two dimensional Hilbert space the performance of quantum feedback control does not depend 
on measurement details, the situation becomes different when we increase the dimension to three. By study- 
ing the three dimensional amplitude damping channel, we have showed there is a class of measurements 
or class of Kraus representations of the amplitude damping channel which perform as good as canonical 
representation. We named it equi-canonical class and analytically found the entanglement fidelity that can 
be attained by performing such measurements. Interestingly, we introduced another class of Kraus repre- 
sentations, the super-canonical class, which leads to an entanglement fidelity higher than the one obtained in 
the equi-canonical class. We analytically found the maximum entanglement fidelity in the super-canonical 
class and the Kraus representations by which this maximum can be attained. By means of numeric tech- 
niques, we discovered that the maximum entanglement fidelity we have found is not only the maximum 
in this class but also the global maximum over all possibilities when considering the most general Kraus 
representations with three Kraus operators. Furthermore our numeric studies shows that the same statement 
is true even if we consider measurement with four outcomes. Indeed, motivated by numerical evidences, 
we feel hke conjecturing that this might be true for any number of Kraus operators (measurement outcomes). 
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Appendix A 

In what follows, we justify the positivity of f2 in (48). Recall that FD{Tcorr) is given by, 



Replacing tr\Do\ and ir-|I?2| from equations (43) and (45) into the above equation and using the following 
identity 




fc=o 



with Di = Ci. Therefore FoiTcorr) becomes, 



FoiTcorr) = ^[{tr\Ci\f + {tr\Do\f + {tr\D2\f]. 



(56) 




(57) 



we obtain 



FniTcorr) = l[tr\\Ci\) + 2p^ + 5(1 -p) + 2,/l^{\j\,/^ + \S\Vk + l)]. (58) 



Using equation (29), we can rewrite the above equation to get equation (47), 
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n 



with 



n = \-f\^/P^h + \5\^/kTl-{2-p). 



To prove that J7 > 0, we first use (57) to rewrite Q as follows. 



n = V + + + + 2MS\Vi9 + h){k + l)-i2-p) 
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= ^(2 - p)2 - 8(1 - p)\jm^ + 2\J\\6\^/p^ + V(l -p) + 16(1 - p)2|7|2|<5|2 - (2 - 
(59) 

Since the following inequality holds, 

2|7ll<5| Vp' + V(l -p) + 16(1 - p)2|7|2|5|2 > 8(1 -p)|7p|<5p, (60) 
we conclude that Q,> 0. 
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